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ABSTRACT 

We evaluate the effect of small scale inhomogeneities on large scale observations 
within the statistics of gravitationally lensed quasars. At this aim, we consider a cos- 
mological model whose large scale properties (dynamics, matter distribution) are the 
same as in Friedmann-Lemaitre models, but whose matter distribution is locally inho- 
mogeneous. We use the well known Dyer-Roder distances to allow a simple analytical 
expression of the optical depth t, and pay particular attention on the different role 
played by the different notions of distances (filled beam angular diameter distance and 
Dyer-Roder distances) when calculating this quantity, following the prescription from 
Ehlers & Schneider for a coherent formalism. We find that the expected number of 
gravitationally lensed quasars is a decreasing function of the dumpiness parameter a. 

Key words: cosmological parameters - cosmology: observations - gravitational lens- 
ing 



1 INTRODUCTION 

\ One of tlie major task in modern cosmology is the precise 
. determination of the parameters which characterize the as- 
sumed cosmological model. In this direct approach (accord- 
ing to Ellis' (1995) terminology) a theoretical description 
' of the space-time is postulated, the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) model, and its parameters are 
determined by fitting the observational data. 

The corner stone of the observational support to the 
FLRW model is the existence and the high isotropy of the 
relic Cosmic Microwave Background Radiation (CMWB). 
Ehlers, Geren & Sachs (1968) showed that if the background 
radiation appears to be exactly isotropic to a given family 
of observers then the space-time is exactly FLRW. There- 
fore, together with the Copernican Principle, we can prove 
the Universe to be FLRW just from our own observations 
of the CMWB. Stoeger, Maartens & Elhs (1995) extended 
this result to the case of an almost isotropic background 
radiation, which implies an almost FLRW space-time. This 
important results is a firmer ground for the assumption of 
the FLRW model to describe the large scale structure of the 
Universe, but it also makes clear that we need to understand 
the departures from a spatially homogeneous model when in- 
terpreting observational data. Indeed, departures form per- 
fect homogeneity change the distance-redshift relation, and 
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this has to be taken into account when fitting the FLRW 
parameters to observations. However, cosmological observa- 
tions are usually fitted just using relationships derived for 
homogeneous models. 

In recent years, several authors have addressed this 
problem in the context of the observations of distant Type la 
supernovae (e.g., Holz 1998, Holz & Wald 1998, Kantowski, 
Kao & Thomas 2000, Sereno et al. 2002, Wang, Holz & Mun- 
shi 2002, Pyne & Birkinshaw 2004). It has been shown that 
the noise due to weak lensing magnification from small scale 
matter inhomogeneities yields large errors on the luminosity 
measurement of high-z supernovae. 

In this paper we aim at investigating the possible sys- 
tematic errors due to neglecting the effects of the local in- 
homogeneities in the distribution of matter when evaluat- 
ing the cosmological constant A from gravitational lenses 
statistics. Statistics of gravitationally lensed multiply im- 
aged quasars has been since a long time considered an useful 
tool to constrain the cosmological parameters, in particular 
the cosmological constant (Turner, Ostrik er fc Gott 1984| 
|Fukugita ot al. 1992 Kochanock 1996), and the properties 
of the lensing galaxies (e.g. Maoz & Rix 1993, Kochanek 
1993). Recently, Chae (2003) has shown that the observed 
gravitational lensing rate in the CLASS radio survey yields 
strong support to a flat cosmological model dominated by 
vacuum energy, with fim — 0.3 and Q.a — 0.7. The precision 
of these results is limited at the moment by the uncertainty 
on the knowledge of the luminosity function of the lens- 
ing galaxy population, their density profile and their evolu- 
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tion since z ~ 1 (e.g., Mao 1991, Chae 2003). For instance, 
Cheng & Krauss (2000) have shown that constraints on cos- 
mological parameter are strongly dependent on the choice 
of galaxies parameters (see also discussion in Kochanek et 
al. (1999)). 

Another major limit is the fact that only few syste- 
matic surveys for multiply imaged quasars have been com- 
pleted up to date fClaeske ns fc Surdej 200T] |; today in fact, 
the statistical uncertainties on Q.a are still dominated by 
the Poisson errors from the small number of gravitationally 
lensed quasars. In the near future, the most promising source 
for new lensed quasars will be wide field surveys (Kuhlen, 
Keeton & Madau 2004) and targeted foUowups of newly dis- 
covered quasars (Morgan 2002, Morgan et al. 2004). For in- 
stance, the Sloan Digital Sky Survey^ will almost double the 
number of known gravitational lenses. Next considerable in- 
crease in the number of gravitational lenses is expected by 
new telescopes like the VLT Survey Telescope^, which will 
allow very wide and deep optical surveys in the Southern 
hemisphere. See also Kuhlen et al. (2004) for a discussion of 
other gravitational lenses surveys to become operational in 
the next future. 

While these large observational projects will consider- 
ably improve the precision of the results and the impor- 
tance of the tool, they also make necessary to consider more 
realistically the detail of the light propagation through the 
observed not homogeneous universe. As pointed out by El- 
lis (1995), small scale inhomogeneities in matter distribu- 
tion have a considerable effect on both observations (Dyer 
& Oattes 1988) and dynamics (Russ et al. 1997) at large 
scale. Moreover, since the lensing effects of the small inhomo- 
geneities on the propagation of light change the angular di- 
ameter distance-redshift relation (Schneider, Ehlers & Falco 
1992), we focus here on this specific problem. Many different 
approaches have been developed to study the gravitational 
lensing in inhomogeneous cosmological models, but the sim- 
plest one from an analytical point of view, and yet efficient, 
is the one proposed by Dyer & Roeder (1972, 1973), in which 
the effect of the local inhomogeneities along the light bun- 
dles are described by the so-called dumpiness parameter a 
(see definition in next section). In the following, we allow the 
dumpiness parameter to be a direction-dependent quantity, 
which is function both of the line of sight to the source and 
its redshift (see, e.g. Wang (1999)). 

The statistical properties of a sample of gravitational 
lenses include the frequency of multiply imaged quasars, the 
distribution of the lenses and source redshifts, of the angu- 
lar separation distribution and of the image multiplicity. In 
this work we will focus on the discussion of the total lens- 
ing probability, leaving a detailed discussion of the other 
statistical properties for a following paper. Moreover, as we 
focus our attention on an effect which is independent from 
our present day knowledge of the galaxy luminosity function 
and the dark matter velocity dispersion, we do not consider 
these aspects in detail. 

The paper is organized as follows. In Sect. 2 we define 
the cosmological model and discuss the relevant distances in 
the study of the propagation of light. In Sect. 3 we calculate 



the gravitational lenses rate, and in Sect. 4 we discuss its 
behavior as a function of a, considering different gravitatio- 
nal lens models. Finally, systematic effects on the estimate 
of the cosmological constant are discussed in Sect. 5, and in 
Sect. 6 we sum up our results. 



2 ROLE OF COSMOLOGICAL DISTANCES 

While cosmological models which are homogeneous at all 
scales are very successful in describing the overall dynam- 
ics and evolution of the Universe, they do not allow a de- 
tailed description of the lensing phenomena. As a fact, all 
the gravitational lensing phenomena (bending of light rays, 
deformation of images, and formation of multiple images) 
are only possible in a clumpy universe, (see, e.g., the discus- 
sion by Krasihsky (1997 1). Therefore, a coherent approach 
needs an inhomogeneous model. 

However, in the statistical analysis of gravitational 
lenses, as well as in analysis of other astronomical obser- 
vations, perfect homogeneity is often assumed (see, e.g., 
discussion in Buchert (2000) about this point). An impor- 
tant reason for this choice is the fact that in homogeneous 
space-times we have simple relations between the proper 
distance and the cosmological distances, i.e. the luminosity 
distance and the angular diameter distance (e.g., Kayser, 
Helbig & Schramm 1997). In inhomogeneous cosmological 
models, these relations are much more complicated, and 
Friedmann-Lemaitre (FL) distances are not generally a good 
approximation to be used in the determination of cosmolog- 
ical parameters from a given set of observational data. 

The relevant distance in gravitational lensing, the an- 
gular diameter distance D, is operationally defined by the 
square root of the ratio of the area dA of a celestial 
body to the solid angle dfi it subtends at the observer 
llSchneider et al. 1992t : 



D = \ —. 



(1) 



In a homogeneous universe, the angular diameter distance 
can be derived from the proper distance Dp using the fol- 
lowing relation 



D{z) = 



l + z' 



(2) 
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This relation does not hold anymore in a clumpy universe. 
The basic reason lies in the fact that the proper distance is 
related to the global geometry of the universe, while the rela- 
tion between the angular diameter distance and the redshift 
is determined mainly by the local matter distribution. 

In this paper we assume a cosmological model which 
is locally inhomogeneous, but homogeneous at very large 
scale, according to some density averaging rule (see, e.g., 
Krasihski (1997)). In other words, we assume that the over- 
all dynamics does not differ from the dynamics of a homoge- 
neous cosmological model. This approximation is well justi- 
fied: by averaging the Friedmann equation, Russ et al. (1997) 
showed that the infiuence of small scale dumpiness on the 
global expansion factor is negligible (so, for instance, the age 
of the Universe can be evaluated using the FLRW relation 
with the Hubble constant, with negligible errors), while the 
distance-redshift relation is significantly affected (so that the 
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measurements of the Hubble constant via measurements of 
standard candles magnitudes have not negligible systematic 
errors). 

In order to describe the degree of inhomogeneity in the 
local distribution of matter, we use a generalized notion of 
the so-called dumpiness parameter (e.g., Dyer & Roeder 
1972, Schneider et al. 1992), introduced by Wang (1999). 
The dumpiness parameter a was first introduced by Dyer 
& Roeder (1972, 1973, hereafter DR), when writing a differ- 
ential equation for the angular diameter distance in locally 
clumpy cosmological models, and it was defined as the frac- 
tion of homogeneously dstributed matter within a given light 
cone. 

Starting from the equations for scalar optics (e.g., Zel- 
dovich 1^564), Kristian & Sachs I|19(i5^ '). DR derived a se- 
cond order difi^erential equation for the diameter of the light 
ray bundle propagating far away from any clumps (i.e. in 
regions where a < 1), assuming a negligible shear. For a 
pressure-less universe with cosmological constant A, the DR 
equation reads: 

(1 + z){n^{i + zY + nK{i + zf + i}A}-^ + 

+ ^an,n{l + zfD = 0, (3) 

with initial conditions D{0) = and -D'(O) — c/Hq. Here 
Qk = 1 — fim — Oa, since we neglect the contribution from 
any relativistic fiuid or radiation. For a = 1 (filled beam 
case) we recover the angular diameter distance, while for 
Q = we have the well-known empty beam approximation. 
For a detailed discussion of the solutions of the equation @ 
within quintessence cosmological models we refer to Sereno 
et al. (2000). For the following. It is useful to introduce the 
dimensionless distance r: 

r(z, fim, f^A, a) = — D(z, //o, rim, riA,a) , (4) 

and the symbol ri for the dimensionless angular diameter 
distance in the filled beam case. 

Note, however that the DR equation ^ is well defined 
for any a > 0, and in its derivation the mass density is never 
required to be uniform. This allows to consider the dumpi- 
ness parameert as a local variable, as done in Wang (1999) to 
descibe the weak lensing magnification of distant standard 
candles. Therefore, in the following we assume the dumpi- 
ness parameter to be a function both of the line of sight and 
the redshift. Given a source at resdhift z in a specific inho- 
mogeneous cosmological model, for any line of sight to the 
observer the dumpiness parameter a is calculated via equa- 
tion where the distance D is given by numerical simula- 
tions. As a consequence, a complete description of the light 
propagation in the clumpy Universe needs the knowledge of 
the probablity distribution function (PDF) for values of the 
dumpiness parameter. We return to this point in Sect. 5 

Let us now consider in more detail the effect of inhomo- 
geneities in observations meant to measure the cosmological 
constant. The DR distance r is a strongly decreasing func- 
tion of a, at fixed redshift (Schneider et al. 1992): there- 
fore, a larger fraction of matter in clumps partially masks 
the effect of a larger cosmological constant when evaluating 
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Figure 1. The Dyer-Roeder distance r(z) for two values of the 
cosmological constant (O^ = 0.0 (continuous curves), 0.7 (dashed 
curves)) in the empty beam case (upper curves) and in the filled 
beam case (lower curves). Space-time is fiat. 



angular diameter and luminosity distances (Fig. 0, since 
there a smaller contribution from the Ricci convergence. On 
the other side, along light beams characterized by a > 1 
(i.e., propagating in overdense regions) agular diameters dis- 
tances are lower than in the filled beam case. For this rea- 
son a large amount of local dumpiness can result in a lower 
value for A when fitting the observational data. In Fig. |5|we 
show the ratio of the DR distance for the empty beam case 
relatively to the filled beam for a few value of the source red- 
shifts. Up to redshift ~ 1, the differences are not important 
in the DR distance itself (although, they maybe not negligi- 
ble for several astronomical observable quantities) . The role 
of the local clumps becomes more and more important at 
higher z. 

Though the assumed cosmological model is not the most 
satisfactory to describe inhomogeneities, since it has not a 
firm theoretical basis in the framework of General Relativity 
(i.e., it is not a solution of the Einstein field equations), it 
allows a simple and efficient description of the light propa- 
gation through a clumpy universe. This model has been dis- 
cussed in detail in Schneider et al. (1992) and Seitz, Schnei- 
der & Ehlers (1994 ). Ehlers & Schneider fl.9g6| have intro- 
duced a self-consistent formalism to study the gravitational 
lensing in a clumpy Universe. In particular, they stressed 
the different roles that the different notions of cosmologi- 
cal distances have in this model. Namely, when we consider 
quantities which are related to the global geometry of the as- 
sumed cosmological model (such as volumes), it is necessary 
to use the FL angular diameter distance. Then, the volume 
element (i.e. the volume of a shell with proper thickness dl) 
is 
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Figure 2. DR distance versus the dumpiness parameter, rela- 
tively to the filled beam case at several value of the redshift. 
= 0.7 in a flat space-time. 



dl/ = 47r ( ) rldl. 



(5) 



In other words, the volume element dV does not depend on 
the local inhomogeneity degree. This is consistent with the 
fact that in the locally inhomogeneous model, on large angu- 
lar scales (i.e., larger than 9 ~ 10", see Linder (1988)), the 
distance-redshift relation is the one computed in the FLRW 
models, for any source redshift. This is also in agreement 
with the results from N-body numerical simulations in Cold 
Dark Matter scenarios (Tomita 1998), where it has been 
shown that the dispersion in values of a along the different 
light paths becomes increasingly larger as the angular is as 
small as a few arcsec. We also note that equation Q can 
be read as the definition of the angular diameter distance in 
homogeneous models (Schneider et al. 1992). 

On the other side, there a simple reason why, when con- 
sidering gravitational lensing phenomena, the DR distance 
has to be used: light deflection modifies the cross area of the 
light bundle and the Ricci focusing term is a linear function 
of the amount of matter within the ray bundle (see, e.g., 
Schneider et al. 1992, sect. 3.4). As a consequence, all the 
distance dependent quantities which play a role in the de- 
scription of the lensing phenomena are functions of the DR 
distance, and, in particular, the strong lensing cross section 
o. 

Let us consider in more detail the cosmological strong 
lensing probability r. This is defined as the probability that 
a hght source at redshift Zs is multiply imaged by a deflector 
at z < Zs- In the expression of r(zs) there are two different 
physical quantities which are functions of combinations of 
distances: the cross section for multiple images a, and the 
volume element dV . 

The cross section a is defined as the area in the lens 



plane for which multiple imaging occurs for sources behind 
it. This quantity depends on the redshift z of the lens and 
a set X of astrophysical parameters which characterize the 
gravitational lens model. For the most generally used mod- 
els (point-like mass distributions and isothermal spheres) 
the cross section depends on a particular combination of 
distances: 



, . J Pis Pi 



(6) 



where Ps-,Pd and Pds are the DR angular diameter dis- 
tances between the observer and the source, the observer and 
the lens and between the lens and the source, respectively. 
We remark that the overall effect of the clumpy distribution 
of matter along the hght rays on the lensing probability is 
not due to any change of the volume element dV , see equa- 
tion lO, but it is only because of the dependence of the 
strong lensing cross section o on the dumpiness parameter. 

In the following we will consider the projection of the 
cross section on the source plane (located at redshift Zb). 



r{zs,a) 



r{z, a) 



(^iz,Zs,x), 



(7) 



where the DR distances are considered. This quantity allows 
a more clear and compact definition of the lensing probabi- 
lity and it is the natural quantity to consider in the assumed 
cosmological model llEhlers fc Schneider 1986t . It is impor- 
tant to note that the quantity a is not in general a function of 
the the distance combination Dds^^d/Ds: as a consequence, 
the point-like mass and isothermal sphere cross sections are 
functions of different distances combinations (see Sect. 4). 
We now calculate the probability of strong lensing phenom- 
ena and evaluate its dependence on a, and then determine 
explicitly the distance functions for these two gravitational 
lens models. 



3 THE STATISTICS OF GRAVITATIONAL 
LENSES 

In this section, we derive the formulae for the statistics of 
gravitational lensing, following mainly the formalism dis- 
cussed in Ehlers & Schneider Ill98fcill . considering in partic- 
ular the proper role of the two types of distances. Let us 
consider a statistical ensemble of cosmological sources at 
redshift Zs and a set of comoving gravitational lenses with 
number density n{z). If we neglect gravitational lenses evo- 
lution and mergings, the comoving number density of lenses 
is conserved: n{z) = no(l + z)'', where no is the local density. 
The number of gravitational lenses in a shell with volume 
dV is then 



N{z) 



n{z)dV 

n{z)A{z)^dz, 



(8) 



where A{z) is the area of the sphere located at redshift z. 

The probability dT{z, Zs,x) that a quasar at z^ is mul- 
tiply imaged by a gravitational lens in the redshift range 
(z, z + dz) is defined as the fraction of the area of the sphere 
at z — Zs (i.e., the fraction of the sky) covered by the gra- 
vitational lenses cross sections cr{z,Zs,x)- This definition 
implicitly assumes that the projected cross sections do not 
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overlap, which is equivalent to state that dT{z,Zs,x) <SC 1. 
The area covered by the projected cross sections of the gra- 
vitational lenses in {z, z + dz) is therefore 



n{z)a{z, Zs, x) A(z)^dz . 

dz 



(9) 



According to the definition, the differential lensing pro- 
bability then reads 



dr(2, Zs) 



n(z)a{z,Zs,x) 



= ■n{z)a{z,x) 



A{z) dl 
A[zs) dz 



dz 



r{zs,a) 


2 


■ ri{z) ' 


r{z, a) 




ri{zs)_ 



^dz. (10) 
dz 



In the particular case a — 1, the distance ratios factor- 
ize, so that we obtain the more common expression 



dr^z, Zs) = n{z)a{z, 



(11) 



i.e., the definition of differential lensing probability in cos- 
mologies which are homogeneous at all scales. Therefore, in 
this case, the relevant distance combination is exactly the 
ratio DdsDd/ Ds which appears in equation Q, and does 
not depend on the particular choice of the gravitational lens 
model. 

Let us now evaluate the explicit expressions for the lens- 
ing probabilities in equation HIH . The quantity dl can be 
written in the following way: 

dl = -cdt = c-^-^, (12) 
rt(z) 1 + 2 

where we considered the past light cone, and H{z) = d/a, 
a = 1/(1-1-2) being the normalized cosmological scale factor. 
In a FLRW cosmological model 

H{Z) = J/o\/!^m(l + 2)3 + (1 - - fiA)(l + 2)2 + f^A .(13) 

Hereafter we focus our attention to flat cosmological 
models, i.e., with SIa + f^m = 1, as these are preferred by 
inflationary scenarios and strongly supported by many re- 
cent observational evidences l |Wang at al. 2000^ . Then, the 
differential lensing probability reads 



dT{z,Zs) = 



"0 7^cr(2,2s,x) 
-no 

(1 + ^)^ 

^n^{i + z)s + nA 



'rizsY 


2 


■ ri{z) ' 


. r{z) . 




ri{zs)_ 



: dz . 



(14) 



It is evident that the properties of the functions dr (and r) 
with respect to a are determined only by the strong lensing 
cross section a{z,Zs). Thereby, it is necessary to consider 
only those functions of the DR distances which enter their 
expressions. In other words, the functions r and a have the 
same qualitative behaviors with respect to the dumpiness 
parameter. Then, our next point is to investigate the be- 
havior of the quantity a{a). 



4 PROPERTIES OF DR DISTANCES 
COMBINATIONS 

In this section we analyze in detail the cross sections of some 
general models of gravitational lenses and the involved com- 
binations of DR distances, in order to evaluate qualitatively 
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Figure 3. The projected cross sections of a point-like gravitatio- 
nal lens for a = 0, 0.5, 0.8, normalized to filled beam case (a = 1), 
versus the lens redshift. Source is at 2^ = 5, and cosmological pa- 
rameters are Ca = 0.7 and Q^n = 0.3 



the dependence of the lensing probability r on the dumpi- 
ness parameter. Asada ( 1998 1 has investigated the analyt- 
ical properties of several combinations of DR distances as 
functions of the dumpiness parameter, and he also deduced 
some consequences on the observable quantities: the deflec- 
tion angle, the time delay, and the lensing probability. Here 
we focus on the study of the quantities which directly enter 
the optical depth. We consider the following models: a point 
mass distribution, singular isothermal sphere, and isother- 
mal sphere with a non zero core radius. In this section we 
only consider the DR distances. 



4.1 Point-like gravitational lens 

A point-like gravitational lens produces two images, what- 
ever the source position is; so, strictly speaking, the given 
definition of strong lensing cross section does not apply here. 
Anyway, it is natural to assume such a cross section (on the 
lens plane) to be the disk with radius equal to Einstein ra- 
dius te. So, the cross section on the source plane reads: 



2 

TrrE 



AtiGM D^Dds 



(15) 



where M is the mass of the lens. This quantity is a decreas- 
ing function of the dumpiness parameter a, as it is shown 
in Fig. 13 Consequently, at smaller value of the dumpiness 
parameter, the lensing probability is higher, when consider- 
ing the lensing cross section for a point-like distribution of 
mass. 

Consider that the distances ratio DdDds/ Ds in the ex- 
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pression of the Einstein radius te is a decreasing function 
of the dumpiness parameter (lAsada 199811 . but it is not 
the right function of DR distances which enters the expres- 
sion of the lensing probability in a DR cosmological model, 
as shown above. As a fact, considering this ratio is not 
coherent with the assumptions on the cosmological model 
IjEhlers fc Sch neider 198611 . and would wrongly lead to pre- 
dict that the lensing probability decreases for decreasing val- 
ues of the dumpiness parameter, as in Asada (il998il . 

4.2 Isothermal spheres 

Let us consider now the gravitational lenses to be isothermal 
spheres, and let us study the dependence of the cross section 
on Q. Many authors have shown that singular isothermal 
spheres (SIS) allow a detailed description of the matter dis- 
tribution in individual gravitational lenses (see, e.g., Rusin 
et al. 2002), and are therefore used in statistical analysis of 
the lensing galaxies population (e.g., Chae 2003). SIS are 
characterized by the surface mass density 

(16) 

where ^ is the position vector on the lens plane, and av the 
line-of-sight velocity dispersion. See, e.g., Hinshaw & Krauss 
H1987[l and Schneider et al. (1992) for a detailed description 
of its lensing properties. The area in the lens plane for mul- 
tiple lensing is 



1671'' 



(17) 



Note that the distance combination which enters equation 
11711 is the same as in the Einstein radius for a point-like 
gravitational lens. But, when we consider the projection on 
the source plane, we get 



16n 



Di 



(18) 



As the angular diameter distance between two points at 

redshifts zi and 22 > 21 is a decreasing function of a 

IjAsada 199811 . so is the cross section cro and, consequently, 

the optical depth r. In Fig. we plot the cross section 

ao{z;a) relatively to the filled beam case, evaluated along 

different line of sghts characterized by different values of the 

dumpiness parameter. 

Let us now consider an isothermal sphere with a 

non zero core radius ^c- The surface mass density is 
2 

S(f) = ^ , i , , which leads to a non constant de- 

flection angle. The cross section on the lens plane is 
lIHinshaw fc Krauss 19871 



5/3-i/?2)-iV^(/3 + 4)«/^ 



/3<i 
/3>i 



(19) 



where /3 is the core radius in units of the natural length 
scale 

>2 DdI>dB 



Co = 4^ — 



(20) 



The projected cross section for multiple imaging is, for f3 < 
1/2, 



where we introduced the quantity C(/3) = a/ao- The pre- 
sence of a core radius lowers the lensing probability, but 
does not change its qualitative properties with respect to 
the dumpiness parameter, see Fig. ^p. Note that the func- 
tion C{P) is not a constant, since (3 is not, with respect to 
the dumpiness parameter. The dimensionless core radius (3 
is a function of the redshifts z,Zs, and the cosmological pa- 
rameters and the dumpiness parameter, via the length unit 
^0: 



1 & 

P{z,Zs,nA,a) = 

47r ai Dd-Ods 



(22) 



Di 



(21) 



It is interesting to note that, despite the fact that the 
cross sections of the considered gravitational lens models 
differ in the DR distances ratio, they are both decreasing 
functions of a. In other words, the monotonic properties 
of the cosmological optical depth with respect to a are ge- 
neral. Finally, it is evident that neglecting a leads, in the 
theoretical predictions, to underestimate the probability of 
strong lensing, and, in the statistical analysis of high redshift 
quasars catalogs, to overestimate the cosmological constant. 
This is discussed in the next section. 



5 ON THE EXPECTED NUMBER OF 
GRAVITATIONAL LENSES 

In this section we evaluate numerically the optical depth in 
a flat cosmology as a function of the cosmological param- 
eter Qa and of the dumpiness parameter a. In particular 
we analyze quantitatively the effect due to neglecting the 
local inhomogeneity in calculating the lensing probability in 
a clumpy universe, along lines of sight with a < 1. For this 
purpose, we consider the gravitational lenses to be singular 
isothermal spheres, because they allow a very simple and 
analytical treatment of the problem, and have been used in 
a variety of studies of the statistical properties of the gra- 
vitational lenses (see Chae (2003) and references therein), 
allowing a direct comparison of the results. Moreover, this 
particular choice does not affect out qualitative results. 

As we are mainly interested in the effects due to a va- 
riation of a, we compare the optical depth t{zs, 0.a, a) with 
To, the numerical value in the case with vanishing cosmo- 
logical constant and light propagation through filled beams, 
so that our discussion is independent of the numerical pa- 
rameter F = 16n^ (ir)^ ("ffe) ^I'^'^l' controls the lensing 
probability. In Fig.|K|we plot the relative optical depth r/ro 
as a function of the source redshift, for six different cosmo- 
logical models, in which = 0.6, 0.8, considering lines of 
sight with Q = 0, 0.5, 1 for each case. As it is well known, 
the strong lensing probability is very sensitive to the value 
of the vacuum energy (e.g., Fukugita et al. 1992). For any 
value of a, the probability that a source at given redshift Zs 
is multiply imaged grows by a factor of ^ 2, if goes from 
~ 0.6 to ~ 0.8. 

In order to disentangle the effect of the clumpy distribu- 
tion of matter from that due to the cosmological constant, 
we consider At/tq, i.e. the relative variation of the lensing 
probability with respect to the case a = 1, as a function of 
f^A along lines of sight characterized by a < 1. We plot this 
quantity, evaluated at different source redshifts, in Fig.|S| As 
shown above, for any Zs and JIa, the lensing probability is 
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Figure 4. Projected cross sections relatively to the filled beam case for a singular isothermal sphere (left-hand panel) and an isothermal 
sphere with core radius = 10 pc (right-hand panel), versus the lens redshift. The source is located at 2a = 5 and cosmological 
parameters are = 0.7, Qm = 0.3. For lens redshift z <^ 1, the two lens models have the same limit. The same qualitative results holds 
for any value of Zs and Qa. 



a decreasing function of the dumpiness parameter. We also 
notice that the effect of the dumpiness parameter increases 
with the redshift of the sources, if a given direction with 
Q < 1 is considered. 

The most important feature to note is that the variation 
becomes rapidly less important at larger values of the cosmo- 
logical constant. The reason for this effect is twofold. First, 
for larger values of the cosmological constant, the influence 
of all other astrophysical and cosmological parameters is ex- 
pected to be less important, since at high Ha, the optical 
depth is very sensible to any small change in the value of the 
cosmological constant. Second, a enters the equation as 
a coefficient of the matter density parameter, determined (in 
the flat cosmological models) by the relation flm = 1 — S7a; so 
it is relatively less important for small values of the density 
parameter. In the most commonly accepted range for the 
cosmological constant (0.6 ^ Qa ^ 0.85), the lensing pro- 
bability increases (relatively to a completely homogeneous 
matter distribution) by a factor of about 7%, 17%, 30% if we 
consider the dumpiness parameter a = 0.75, 0.5, respec- 
tively. 

As stated above, in order to have a coherent description 
of the small scale dumpiness we need to consider the PDF 
for a, since this a direction-dependent quantity. The PDF 
depends on the background FLRW cosmological model, and 
can be derived from the PDF of the lensing magniflcation 
^J.{z), calculated in numerical simulations (e.g., Wambsganss 
et al. 1997, Holz & Wald 1998, Tomita 1998, MortseU 2002), 
via the relation 



and the DR equation. As shown by Wang (1999), in general 
the a PDF is peaked at a < 1 for any z, but it tends to be 
more symmetric and shows smaller scatter around the peak 



value as the redshift increases (see, e.g., Fig. 3 in Wang 
1999). In Fig. |7| we plot, for a flat spacetime with flm = 
0.4, both the most likely value and the average values, as 
calculated via the approximate analytic expressions of the 
PDF (q) given in Wang (1999). 

The mean value of a is 1 at any redshift in all cosmo- 
logical models (if they are FLRW on average as assumed 
here). This is indeed the same proprierty which is known to 
hold for the PDF of the lensing magnification fJ,{z), whose 
basic motivation is the flux conservation (Weinberg 1976). 
Therefore we expect that the effect of inhomogeneoties on 
the lensing cross sections is reduced in large ensemble of 
gravitational lenses, or very distant sources. 

For a gravitational lens at z ~ 1, the most likely value is 
Q ~ 0.85; this translates in an underestimates of the lensing 
cross section for a SIS (using the fiUe beam distances and 
considering a source at Zs > 4) of a factor ~ 1.15, see Fig.|l| 
These systematic errors decrease for more distant gravita- 
tional lenses or sources, but since the most likely value is 
always lower than 1, and in a relatively small sample of gra- 
vitational lenses the mean of a nonsymmetric probability 
distribution is not likely the best estimator, such errors may 
become not neglibible when evaluating the cosmological pa- 
rameters. 

While a detailed statistical analysis including an accu- 
rate determination for the PDF for any cosmological model 
is beyond the scope of the present work, it is anyway impor- 
tant to estimate now an upper limit to the possible errors on 
the predicted number of gravitational lenses. We have there- 
fore considered the list of 1163 luminous quasars^, which 
have been observed in the following optical surveys; CFHT 

^ This catalog is available at the web address 
|http : / /vela ■ astro ■ ulg ■ ac ■ be/themes/extragal/gravlens/bibdat 
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Figure 5. The optical deptli as a function of the redshift for 
differ values of a and SI a relatively to tq, the optical depth in 
the case Oa = and a = 1. Continuous and dashed curves are 
for Qa = 0.6, 0.8 respectively. The dumpiness parameter is a = 
0, 0.5, 1.0, from the above curve to the lower one for each value of 
SIa- Gravitational lenses are SIS and space-time is flat. 

(Crampton, McClure & Fletcher 1992), CFHT (Yee, Filip- 
penko & Tang 1993), HST (Maoz et al. 1993), NOT (Jaunsen 
et al. 1995). This catalog contains 7 confirmed gravitational 
lenses, and its redshift distribution is plotted in Fig. |S| the 
peak is at redshift 2; ~ 2, and only a small fraction of sources 



Figure 7. The peak (continuos line) and the mean value (dashed 
- dotted line) of the dumpiness parameter in flat cosmological 
model with £7^ = 0.4, calculated via the analytic approximation 
given in Wang (1999). 

are beyond 2 = 3. We have calculated the relative variation 
in the expected number of multiply imaged quasars consid- 
ering different values for the dumpiness parameter, corre- 
sponding to peak values of its PDF for z ~ 0.5 and 3.0. 
The lensing galaxies are being modelled as SIS (note that 
including a small core makes the dumpiness effect slighlty 
larger). In Fig. |^ we plot the relative variation of the ex- 
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redshift 



Figure 8. Redshift distribution of the 1163 luminous quasars in 
the catalog used to calculate the effect of a dumpiness parameter 
a 7^ 1. See text for details. 

pected number of lenses as a function of Qa for the peak 
values a = 0.75, 0.95. this figure shows the upper limit for 
the systematic errors that can be found for the predicted 
numner fo gravitational lenses when adopting the simple 
filled beam hypothesis; given the property a = 1, this could 
be further reduced in upcoming larger surveys. Note that, 
since we are considering a fiat spacetime and the effect of the 
local inhomogeneities increases with Qui, the variation of the 
expected number A'' of multiply imaged quasars is a decreas- 
ing function of the cosmological constant. However, only for 
a very large value of the cosmological constant {Qa > 0.8), 
such a variation becomes rapidly small. This makes clear 
the point that if in the evaluation of r we use angular di- 
ameter distances for a perfectly homogeneous cosmological 
model, we can underestimate the lensing probability and, 
consequently, overestimate Qa- 



6 CONCLUSIONS 

In this work we have investigated whether the local depar- 
tures from a completly homogeneous cosmological model can 
have observable effects in the statistical study of high-z gra- 
vitational lenses. Following the work by Ehlers & Schneider 
1198bl . we derived the expressions for the cosmological opti- 
cal depth in the framework of a cosmological model which is 
FLRW on very large scales (i.e., whose overall dynamics is 
very well described by FLRW models) and whose matter dis- 
tribution is locally inhomogeneous. The direction-dependent 
dumpiness parameter a quantifies the fraction of matter in 
compact objects along a given line of sight, and its peak and 
mean values (as a function fo the source redshift) are cal- 
culated via the analytical approximation of the PDF given 
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Figure 9. Relative variation of the expected number of strong 
lensed quasars in the given sample of high luminosity quasars with 
respect to the filled beam C3.S6j ciS Si function of the cosmological 
constant. The dumpiness parameter is a = 0.75 (dashed curve) 
and a = 0.95 (continous curve). Lenses are modelled as SIS. 



in Wang (1999). We have paid particular attention to dis- 
entangle the different role played by the different notions of 
distance in the definition of the optical depth T{a): the small 
scale inhomogeneities along the line of sight do not change 
the volume element dV, but the strong lensing cross section, 
see equation ||KJ and Fig.|l] 

Up to redshift z ~ 3, the most probable value of the 
dumpiness parameter is very different from the average 
value (which is constrained to be 1 in a model homoge- 
neous on a large scale), see Fig. |7| and the effect may be 
important in statistical analysis of relatively small sets of 
gravitationally lensed sources. 

Asada (1998) presented a similar calculation to the 
one done here in Sect. 4. He discussed the influence of the 
dumpiness parameter on several gravitational lensing ob- 
servable quantities, showing that in a clumpy universe de- 
flection angles are smaller and time delays are longer than in 
an homogeneous universe, given the same lens-source con- 
figuration. However, in contraddiction with Asada's result, 
we have found that the gravitational lenses rate is a decreas- 
ing function of the dumpiness parameter. We have shown 
that in the empty and filled beam cases different angular 
distances ratios enter the optical depth expression (Sect. 3), 
and this leads to an higher number of expected of gravita- 
tional lenses when light beams with a < 1 are considered. 

While a detailed statistical analysis including the effect 
of small scale inhomogeneities on the determination of the 
cosmological constant is beyond the scope of the present 
work, we can already draw some important conclusions and 
compare our findings with previous works. Since the a mean 
value tends to be 1 when a sufficiently large number of dif- 
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ferent lines of sight is considered, the effect described here 
tends to be less important in large surveys for gravitation- 
ally lensed sources. Moreover, since the peak value also tends 
to 1 for z > 5, statisical study of very high-redshift sources 
is less affected by the local dumpiness along the different 
lines of sight. On small set of lensed sources, since the peak 
value of the dumpiness parameter is always lower than 1, 
using the filled beam angular diameter distances leads to an 
overestimate of the number of the expect gravitational lenses 
(Fig. 9), and, consequently, to overestimate the cosmological 
cosntant. 

It is interesting to note that the same qualitative re- 
sult (on the determination of the cosmological constant) has 
been found in previous works in which the authors have 
taken into accounts the effect on the distance-redshift rela- 
tion playd by the matter distribution inhomogeneities. Kan- 
towski (1998) described the effects of inhomogeneities on the 
determination of the cosmological constant and i},n using a 
"Swiss cheese" model to derive analytic expressions for the 
distance-redshift relation. He has shown that, when analyz- 
ing high-z standard candles, assuming a completely homo- 
geneous matter distribution leads to overestimating the cos- 
mological constant if the filled beam hypotheisis is always 
used (see, e.g.. Fig. 8 in Kantowski (1998)). 

In a recent work, Barris et al. (2004) analyze a set of 
194 Type la SNe at z > 0.01 (including the 23 discovered 
in the IfA Deep Survey in the range z = 0.34 — 1.03) us- 
ing both the filled beam and the empty beam assumption 
for all the lines of sight to calculate luminosity distances^ . 
Although the presence of small scale inhomogeneities is far 
from eliminating the need of a vacuum energy contribution 
to explain the new data, they do change the final estimate 
of the cosmological parameters flm and JIa: indeed, when 
using the angular-diameter distances with a = 0, the confi- 
dence contours for the comsological parameters are shifted 
to lower values of Ha and higher values of flm 

The same effect has been found by authors which inves- 
tigated the effect of large scale inhomogeneity in the matter 
distribution. Tomita (2001) considered a local void on scales 
of 200-300 Mpc around our Galaxy to interpret the high- 2: 
SN data, effectively describing it with with a < 1 in the DR 
distances. He finds that the data (available at the time) coule 
be well fitted with Qa ~ 0.4 (value estimated in the over- 
dense outer region). It is certainly interesting to test these 
findings with the more recent and larger set of data available 
at the moment. In conclusion, a precise determination of the 
cosmological parameters using both the Hubble diagram of 
the high-2 Type la SN and the statistical properties of gra- 
vitational lenses requires an accurate determination of the 
effect of the local dumpiness in the matter distribution on 
the light propagation. 



^ Note, however, that using the empty beam approximation for 
all the lines of sight is not a coherent description, since in an 
inhomogeneous Universe not all the light beams can be devoided 
of matter, and more realtic distribution for a has to be chosen to 
describe the dumpiness effects using the DR, distances. 
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